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Abstract 

We introduce an event-by-event pQCD -r saturation -l- hydro (”EKRT”) framework for high-energy heavy-ion col¬ 
lisions, where we compute the produced fluctuating QCD-matter energy densities from next-to-leading order (NLO) 
perturbative QCD (pQCD) using saturation to control soft particle production, and describe the space-time evolution 
of the QCD matter with viscous hydrodynamics, event by event (EbyE). We compare the computed centrality depen¬ 
dence of hadronic multiplicities, pr spectra and flow coefficients Vn against LHC and RHIC data. We compare also 
the computed EbyE probability distributions of relative fluctuations of Vn, as well as correlations of 2 and 3 event- 
plane angles, with LHC data. Our systematic multi-energy and -observable analysis not only tests the initial state 
calculation and applicability of hydrodynamics, but also makes it possible to constrain the temperature dependence 
of the shear viscosity-to-entropy ratio, ri/s{T), of QCD matter in its different phases. Remarkably, we can describe 
all these different flow observables and correlations consistently with rj/s{T) that is independent of the collision 
energy. 


1 Introduction 

The basic idea in the new NLO-improved perturbative QCD + saturation + viscous hydro event-by-event EKRT frame¬ 
work ID presented in this talk, is as follows: Using collinear factorization and NLO pQCD, we first compute the 
production of transverse energy {Ep) carried by partons of pp few GeV (minijets) into the central region in A+A 
collisions. The conjecture of gluon saturation, implemented as in the NLO-extension 121 of the original EKRT paper 
13, then allows us to compute QCD-matter initial conditions, EbyE, for viscous hydrodynamics with which we de¬ 
scribe the space-time evolution of the produced QCD matter. Comparing with a multitude of experimental RHIC and 
LHC data for low-p^ observables, we aim at (i) pinning down the rj/s{T) ratio, (ii) testing our initial state calculation, 
and (Hi) studying the applicability region of viscous hydrodynamics. 


2 Minijet Et production from NLO pQCD 


Using collinear factorization, NLO minijet Ep production per transverse area into a rapidity interval Ay in high-energy 
A+A collisions of impact parameter b can be computed as 


^ = Ta{v + h/2)TA{v - h/2)a{Ep), 


( 1 ) 


where the collision geometry is given by the the standard nuclear thickness functions Ta- The minijet part a, which 
in NLO 15] m [D can be written as 


r{Ep) = X U? 


d|PS1.5j5^& 


( 2 ) 


contains phase-space integrated 2^2 and 2 —3 differential ) partonic cross sections dcr^“*^"/d[PS]^, obtained 
in terms of UV-renormalized squared matrix elements l6| |71, and CTEQ6M IS] parton distribution functions with 
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EPS09s spatially dependent nuclear effects 0. The measurement functions, 


> 2po^ X S 6{yi G Ay)pTi > , (3) 

contain step functions 9 and are analogous to those in jet production Go). They define the 1) minijet Et as the sum of 
the minijet py’s in Aj/; 2) pT cut-off scale po ^ Aqcd above which we do the minijet calculation, and 3) minimum 
Et that we might require in the interval Ay. Such ensure a well-defined NLO calculation safe from collinear and 
infra-red singularities. With these ingredients, the minijet Et calculation is field-theoretically rigorous for hxed pq 
and P G [0,1]. As the Et here is not a direct observable, we must leave the parameter /3 to be hxed from experimental 
data. 


Sn = i^0{y,eAy)pn 


3 Saturation in NLO minijet production 

The parameter po controlling minijet production is here hxed based on a gluon saturation conjecture. As argued in 
El, the Et production is expected to cease (saturate) when 3 —?> 2 and higher order processes become of the same 
magnitude as the basic 2 —2 ones, gi^(2 —?> 2) ~ g7^(3 —>■ 2). This leads to a saturation criterion analogous to 
that in the original EKRT-model ||3l but now for an IR/CL-safe Et, 

VsJvtv.A, Ay,r,b;/3) = ^^plAy, (4) 

d^r TT 

from which po = Psatiy/ snn, A, Ay, r, b; /3, ATgat) is solved locally in r and for hxed /?, ATgat- The ’’packing fac¬ 
tor”, proportionality constant K^at, is to be hxed from the data. The key observation EH [Ill enabling the EbyE 
framework here, is that psat scales essentially as [TaTa]^, so that we can parametrize Pgat(TATA) vs. /3, Kgat- The 
parametrization, intended also for public use, can be found in El- The minimum pgat we allow here is 1 GeV. 


4 NLO EKRT EbyE framework 


Eor the EbyE set-up, we hrst sample the nucleon positions in the colliding nuclei from the standard two-parameter 
Woods-Saxon density. Around each nucleon, we then set a gluon transverse density, T„(r) = (27r(T^)“^ exp(—r^/(2cr^)), 
whose width parameter, a = 0.43 fm is obtained from the exclusive electroproduction data of J/ij) at HERA/ZEUS 
E3]. After this the thickness functions Ta{x, y) of the nuclei can be computed, and pgat (TaT^) obtained for hxed 
P, ATgat on the basis of the NLO pQCD-tsaturation calculation described above. Erom p^at, we then extract the local 
formation time, rs(r) = l/psat(r), and energy density 


e(r,rg(r)) 


dET 1 
d^r rg(r)Ay 


Aigat 

TT 


[Psat(r)]^ 


(5) 


of the minijet plasma. 

”Pre-thermal” evolution from Ts(r) to tq = 1/p™” = 0.2 fm is here done simply with 1 D Bjorken hydro (also 
free streaming would work EH). At the edges of the collision system, below p™" = 1 GeV, we connect the e-prohle 
smoothly to a binary prohle (see El). 

With such EbyE initial conditions, we then run 2 h- 1 D viscous hydrodynamics, EbyE. Our hydro setup is the 
one used in El ES Uhl [HI, i.e. 2nd-order dissipative relativistic hydro with transient huid-dynamics equation of 
motion for the shear-stress tensor EmSKH (see again El)- We neglect heat conductivity and bulk viscosity 
and study the effects of shear viscosity only. We explore the possible temperature dependence of r]/s{T) with the 
parametrizations shown in Eig. 

As the QCD-matter equation of state, we employ the s95p-PCE-vl parametrization EOl with chemical decoupling 
at a rather high temperature of 175 MeV, and kinetic freeze-out at 100 MeV. Resonance decays after the freeze-out are 
included. The initial and transverse flow vt{tq) are set to zero. In the non-equilibrium particle distributions 

on the freeze-out surface, we assume, as usual, that the relative deviations from the equilibrium distributions are 
proportional to p^pyTr^'^ for each particle species. 
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Figure 1: The tested temperature dependences of rj/s. From (T). 


5 Comparison with LHC and RHIC data 

Let us then consider selected results from the extensive multi-energy and multi-observable analysis of Ref. IT]. Figure 
1^ shows the centrality dependence of the computed charged hadron multiplicity dNch/drj and their flow coefficients 
u„{2} from the 2-particle cumulants at the LHC and t;2,3{2}, t;4{3} at RHIC, and comparison with the data. For 
details, definitions of these quantities, and pT spectra, see ID. For each r]/s{T) parametrization, and setting f3 = 0.8 
CD, we exploit only the centralmost LHC multiplicity datapoint to fix the normalization with iCgat- The good agree¬ 
ment with the data indicates that our initial states are on the average working well. All the rj/s{T) parametrizations 
of Fig. [^reproduce the measured r;„{2}’s at the LHC while the RHIC data favors a small p/siT) in the hadron gas 
phase. 

Figure shows the probability distributions of the charged hadron 5v2 = {v 2 — {V 2 ) ev) /{^ 2 } ev and the initial 
spatial eccentricities Se 2 = {€2 — {^ 2 )ev)/{^ 2 )ev at the LHC in different centrality classes. Our EbyE framework 
catches remarkably well the centrality systematics of the measured P{6vn), and also proves that since i5e2 and Sv 2 
are nonlinearly correlated in non-central collisions, hydrodynamics (collectivity) is indeed required to reproduce the 
measured P{Svn)- Eurthermore, the P{Svn) offer direct constraints for the initial state in that they are insensitive to 

vis- 

Eigurej^shows examples of correlations of two and three event-plane angles, (cos(/cifE'i + ■ ■ • -|-nfc„'k„))gp at the 
LHC. Especially since the P{Svn) now constrain our initial state calculation independently of p/s, these correlations 
give vital further constraints for the viscosity and for the validity of the EbyE viscous framework. Remarkably, the 
same two p/s{T) parametrizations that explain the RHIC t;„s (black and blue in Eig.[^ work best also at the LHC! 
Importantly, as the rightmost panel shows, for these cases also the 6f corrections remain conveniently small from 
central to semicentral collisions, suggesting that the obtained hydrodynamic results appear trustworthy up to 40-50% 
centrality classes. Eor more data and discussion, consult again Ref. ID- 

6 Conclusions and outlook 

Performing a simultaneous LHC and RHIC multi-observable analysis, we have shown that the new NLO-improved 
EbyE EKRT framework explains remarkably consistently the LHC and RHIC bulk observables in URHIC. The frame¬ 
work has clear predictive power in cms-energy and centrality, and it is a promising tool for getting a controlled 
estimate of the T-dependence of the QCD matter shear viscosity, now that enough orthogonal data constraints start to 
be available. Our results favor a QCD-matter p/s{T) which is between a constant value 0.2 and a modestly decreas- 
ing(increasing) in the hadron gas(QGP) phase. Similar magnitudes have been obtained also for constant p/s by other 
groups EollSD. Next, we need a genuine global analysis to extract true statistical error limits to p/s{T). Also the 
effects of bulk viscosity should be studied. Work for the EKRT EbyE framework predictions for the forthcoming 
5 TeV PbH-Pb run at the LHC is in progress. 
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Figure 2: Upper panels: Centrality dependence of charged hadron multiplicities at the LHC (a) and RHIC (b). Lower panels: 
Centrality dependence of the 2-hadron cumulant flow coefficients Vn{2} at the LHC (a), and V2{2}, V3{2}, V4{3} from the charged 
hadron 2- and 3-particle cumulants at RHIC (b). Experimental multiplicity data are from flli (a) and (22l|^ (b), and Vn data from 
M (a) and |25l|26l|27l (b). From d. 



Figure 3; Probability distributions of the charged hadron Sv 2 , and of the initial S £2 (and 5ei,2, see (T)) in three different centrality 
classes in ^sjvjv = 2.76 TeV Pb-l-Pb collisions at the LHC. Figures from |[T|> and experimental data from ATLAS 1281 . 
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Figure 4: Examples of correlations of two {left) and three (middle) event-plane angles for charged particles at the LHC, compared 
with the ATLAS data 1291 . Color code as in Fig.[^ Right. As the left panel but computed for p/s =0.1 (black), 0.2 (blue), 0.3 
(red) and parami (orange) with (solid) and without (dashed) Sf corrections. From JT(- 
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